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ABSTRACT 

Shocks and blastwaves are expected to be driven driven into the intracluster medium filling galaxy groups and 
clusters by powerful outbursts of active galactic nuclei or quasars in the member galaxies; the first footprints 
of shock fronts have been tentatively traced out with X-ray imaging. We show how overpressures in the blasts 
behind the shock can prove the case and also provide specific marks of the nuclear activity: its strength, its 
current stage, and the nature of its prevailing output. We propose to detect these marks with the aimed pressure 
probe constituted by the resolved Sunyaev-Zel'dovich effect. We compute and discuss the outcomes to be 
expected in nearby and distant sources at different stages of their activity. 

Subject headings: cosmic microwave background - galaxies: clusters - quasars: general - shock waves 
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1. INTRODUCTION 

Density jumps have been recently pinpointed by X-ray 
imaging of the hot intracluster medium (ICM) that pervades 
galaxy groups and clusters with average densities around 
n ~ 10 cm" 3 and temperatures T in the keV range. These 
jumps have been interpreted in terms of shock fronts propa- 
gating into the ICM out to radial distances r « 0.2 Mpc, with 
Mach numbers around A4 « 1.5, and involving energies up to 
AE «3x 10 61 ergs (Mazzotta et al. 2004; Mc Namara et al. 
2005; Forman et al. 2005; Nulsen et al. 2005a, b). 

Shocks over large scales with such intermediate strengths 
were specifically expected by Cavaliere, Lapi & Menci (2002) 
as marks of the energy being fed back into the ICM by active 
galactic nuclei (AGNs) when they flare up in member galax- 
ies of a group or cluster. Such events occur when a central 
supermassive black hole (BH) accretes an additional mass 
M, ~ 10 9 M Q ; with standard efficiency rj ~ 10 _1 for mass- 
energy conversion, this yields over times Af ~ 10 8 yr energies 
of order 2 x 10 62 (M./10 9 M Q ) ergs. 

If these outputs couple at levels / ~ a few percents to the 
surrounding ICM, they constitute impulsive, considerable ad- 
ditions AE m 10 6I (//5%)(M./10 9 M Q ) ergs to its binding 
energy. In fact, the latter comes to E « GMm/Ar « 3 x 10 61 
ergs in the central 0.2 Mpcs of a cluster, that encompass a 
dark matter (DM) mass M~5x 10 I3 M o and an ICM mass 
fraction m/M » 0.15. With such appreciable ratios AE/E, 
we expect in the ICM a large-scale blastwave bounded by a 
leading shock that starts from the host galaxy and moves into 
the surrounding ICM out to several 10 _1 Mpc. 

In fact, the shock Mach numbers are provided in the simple 
form M w (1 + AE/E) 1 ! 1 by the hydrodynamics of the ICM 
(see Lapi, Cavaliere & Menci 2005), a good electron-proton 
plasma which in thermal equilibrium (see § 4) constitutes 
a single "monoatomic" fluid with pressure P w 2nkT. For 
BH masses bounded by M. < 5 x 1O 9 M (Ferrarese 2002; 
Tremaine et al. 2002) relative energy inputs AE/E < 1 ob- 
tain and yield just M. w 1 .5; they also yield standard Rankine- 
Hugoniot jumps of the post- to the pre-shock density «2/wi = 
4M 2 /(M 2 + 3) « 1.7, consistent with the X-ray analyses. 
What other marks will establish such shocks and blasts? 
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One is constituted by the temperature. This rises sharply 
across a shock but then falls down in the blast, if nothing 
else by adiabatic cooling. In any case, resolved spectroscopic 
measurements of the electron T require many X-ray photons, 
more than currently available from distant clusters or groups. 

Pressure provides another, independent mark; the electron 
pressure is directly sensed with the SZ effect (Sunyaev & 
Zel'dovich 1972). The pressure at the shock is to jump up 
from the unperturbed level P\ by the factor P2/P1 = (5 M 2 - 
l)/4, which marks a shock from a cold front; moreover, P 
must retain sufficiently high levels throughout the blast as to 
propel forward the ICM it sweeps up. 

We will see that the radial pressure run P(r) actually rises 
from the leading shock to an inner "cavity", as long as the 
blast is driven on by a central AGN. Power must be transmit- 
ted from it to the surrounding ICM blast by means of an inter- 
vening medium. This may be constituted either by relativistic 
particles filling up a radiovolume energized by jets (Scheuer 
1974; Heinz, Reynolds, & Begelman 1998); or by another 
and hotter plasma heated up by the impact of radiation-driven 
superwinds (see Lamers & Cassinelli 1999). 

Thus SZ pressure probing will have direct implications for 
the kind and the time-cycle of the AGN outputs, in particular 
of the radio-loud components. 

2. COMPUTING OVERPRESSURES IN BLASTWAVES 

In predicting pressure distributions, a divide is set by the 
source time scale Af compared to the blast crossing time 
Rs/Rs ~ 2 x 10 8 M~ l (# s /0.2Mpc) yr; heretf, is the shockra- 
dial position, and R s = M. cs its velocity in terms of the sound 
speed cs = (5P/3nm p ) 1 ' 2 . 

The impact on the ICM of a short lived AGN source may be 
modeled as an instantaneous central explosion launching the 
classic self-similar blasts of Sedov (1959) and Parker (1963); 
these propagate freely at high if decreasing Mach numbers M. 
(which make Pi and gravity irrelevant) into an initial density 
run n(r). When the latter is provided in the form n(r) oc r~ 2 by 
the standard atmosphere in isothermal equilibrium, the pres- 
sure P in the blast declines toward the center. Steeper density 
runs yield an even stronger decline, whilst flatter ones still 
imply P <P 2 . 

But in some sources the observations recalled in § 1 show 
proximity of the shocks to the edge of the radiovolumes. This 
indicates blasts being currently driven by the relativistic parti- 
cles associated with the radiosource, and calls for considering 
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AGN outputs sustained during the propagation. A similar in- 
dication is provided by the evidence at redshifts z > 0.3 of su- 
perwinds driven by shining quasars (see Stockton et al. 2006). 
The moderate values observed for M require including the ef- 
fects of gravity and finite initial pressure Pi . 

These conditions are aptly modeled with the self-similar 
blasts derived and discussed by Lapi et al. (2005). A sig- 
nificant testbed is again provided by the isothermal distribu- 
tion of both the DM and the unperturbed ICM in the gravi- 
tational field of the former; their cumulative masses out to r 
scale like m{< r) oc M{< r) oc r, and the ICM pressure like 
P\{r) oc ni(r) oc r~ 2 . 

The overall parameters of the perturbing blast may be de- 
rived from the simple "shell approximation", long known to 
be effective and precise (see Cavaliere & Messina 1976; Os- 
triker & McKee 1988). This treats the blast as a shell con- 
taining the mass m{< R s ) swept up out to the radius r = R s , 
and propelled outward by the volume-averaged pressure (P) 
against, we add, the upstream pressure Pi and the DM gravity 
GM{< R s )m{< R S )/R 2 S . It leads to the momentum equation 

j t lm« R s )v 2 ] = 4nR 2 [(P) -Pi] - ° M{ < Rs) m{< R s ) (1) 

in terms of the postshock speed v 2 = 3(M 2 - \)R S /4M 2 . 
Here the pressure terms scale following PR 2 = const, as does 
the gravitational term. Self-similar solutions require also 
the l.h.s. to remain constant as the blast moves out; since 
m{< P.,)v2 oc R S R S - R s cj/R s holds (with c s = const), the 
only consistent solution must satisfy R s oc c S - Thus the shock 
moves outward following R s = M. c s t . 

On the other hand, the gravitational energy in the DM po- 
tential well scales like GM{< R s )m{< R s )/R s oc R s , and so 
do all energies including the binding E; this implies for the 
energy injection AE(t) oc t to hold, that is, the power L(t) = 
const is sustained over a crossing time. Then the blast runs 
unattenuated with AE(t)/E(< R s [t]) andM w (1 + AE/E) 1 ! 2 
independent of time and position. 

The actual existence of the self-similar solutions is proven 
from the full hydro equations. By self-similarity these reduce 
to a set of three ordinary differential equations in the normal- 
ized variable r/R s , that can be integrated inward of the shock 
where Rankine-Hugoniot boundary conditions hold, see Lapi 
et al. (2005). Here we focus on the shape of P{r/R s ) that is 
plotted in detail in Fig. 1 {left panel). 

The point to stress is that P(r/R s ) rises from the shock to the 
blast inner boundary, namely, the "piston". This is the contact 
discontinuity separating the ICM from the inner medium in 
the cavity; the mass m{< R p ) of swept-up gas piles up there, 
causing the density to diverge mildly whereas the temperature 
T{r) vanishes. In the blast the two combine to yield a rising 
P{r/R s ); at the piston this attains a value P p exceeding P2 by 
factors 3-4 (see Fig. 1), larger at lower AL 

The rising trend is best understood in the shell approxima- 
tion on computing the coefficients of Eq. (1), given the scal- 
ings above; simple algebra and some labor yield 
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The first additive term in the middle expression embodies the 
contribution to (P) from the shell momentum; the second, 
those from the initial pressure (contributing 1/3) plus grav- 
itation (2/3). The final expression shows that (P)/p2 and 



hence the monotonic P(r)/P 2 go to unity for M. S> 1, the re- 
sult known from Sedov blasts launched into a density gradi- 
ent with negligible gravity and initial pressure. On the other 
hand, (P) and P(r)/P 2 increase as M —> 1; this concurs with 
the lowering jump P2/P1 to yield an overall value P p which 
decreases slowly with AL as shown by Fig. 1. 

In sum, a driving source excavates in the ICM a cavity with 
very low inner X-ray emission; this is bounded by the blast 
(from a rim at the piston to the leading shock) in a cocoon- 
like topology. Enhanced values Pp/Pi relate mainly to gravity, 
important to the ICM as long as all velocities do not exceed 
cs, or the DM velocity dispersion; so at lower Af a relatively 
stronger push is required to propel the cocoon on. 

From the energy standpoint, high values of P p may be 
viewed as necessary for soaking up - in the form of work 
4nR 2 P p R s done at the piston - the power fraction transferred 
from the source via the inner medium filling the cavity. The 
fraction is a « 1/2 when the medium is constituted by rela- 
tivistic particles, see Cavaliere & Messina (1976), and a = 1 
in the case of a thermalized superwind. 

3. MEASURING OVERPRESSURES WITH THE SZ EFFECT 

The enhanced pressures in the blast can be directly probed 
with the SZ effect. This arises when photons of the cosmic 
microwave background (CMB) crossing a cluster are Comp- 
ton upscattered by the hot ICM electrons; then the black body 
spectrum of the CMB is tilted slightly toward higher energies. 

The resulting intensity change Al of the CMB radiation at 
a normalized frequency x = hv/kTcus is given in the thermal, 
non-relativistic case (Rephaeli 1995) by 



a , o (£?cmb) 3 , . 
A/ = 2 g(x)y. 
{hey 



(3) 



The effect strength is set by the Comptonization parameter 

^max 



y(w) = 2 



m e c 2 



d£p(r) , 



(4) 



just proportional to the electron pressure p = nkT integrated 
along the line of sight (l.o.s.) at a projected distance w 
from the cluster center; £ max = (R 2 + w 2 yl 2 applies if the ICM 
boundary is set at the virial radius R. 

The spectral shape is encoded in the factor g(x). This has 
a crossover point at v k, 220 GHz; it is positive for larger v 
with a peak at v m 370 GHz, and negative for lower v with a 
minimum at v s=a 130 GHz. The approximation g(x) ~ -2x 2 
applies at the Rayleigh-Jeans end of the CMB spectrum. 

As Al oc y oc P holds, we expect that in an AGN-driven 
blast the SZ signal is enhanced from shock to piston, relative 
to the equilibrium value A/ eq on the l.o.s. grazing the shock. 
In Fig. 1 {right panels) we show the outcome of our computa- 
tions based on Eq. (3) and on P(r/R s ) given on the left. 

We describe our results beginning with v < 220 GHz. The 
rise of A//A/ eq shown in the top right panel starts up at the 
shock position w = R s with a high derivative arising from the 
contribution of the pressure jump P2/P1 to the l.o.s. integral. 
The subsequent rise toward the piston reflects the further pres- 
sure increase into the blast. 

Inward of the piston A//A/ eq rises still further when the 
inner medium is constituted by a very hot thermal plasma with 
constant pressure P p {dashed lines); this is closely the case 
with a thermalized superwind, where c s is high owing to the 
high pressure P p ^> Pi and the low density n <n\. 

On the contrary, Al/Al eq declines for w < R p {solid lines) 
when the cavity within R p is filled with relativistic electrons 
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FIG. 1. — Left panel: overpressure in an AGN-driven blastwave; top right panel: enhanced SZ signal (relative to A/ e q on the l.o.s. grazing the shock) at an 
instrumental resolution of 1" in a cluster at z = 0.1; bottom right panel: same for a resolution of 10". Thick lines are for Mach number M = 1.5, when the 
piston location is at R p xs 0.62R,; and thin lines are for M = 2, when the piston is at R p xs 0.72R,. Solid lines are for a blast continuously driven by relativistic 
particles; dashed lines are for a blast continuously driven by a hot inner plasma; dotted lines are for a standard Sedov-Taylor blastwave launched by a sudden 
energy release. 



(with cs ~ c/y/5). In fact, these provide the pressure P p w 8 Pi 
but contribute little to the overall SZ signal on two grounds. 
First, relativistic pressures p K \ = P p imply low electron densi- 
ties even with minimal Lorentz factors of a few (see Pfrom- 
mer, EnBlin & Sarazin 2005; Colafrancesco 2005). Second, 
relativistic electrons are generally inefficient contributors at 
a given v since the tilt they cause in the CMB spectrum is 
stretched toward high frequencies. 

For comparison, we also show the thermal SZ effect pro- 
duced by free Sedov-Taylor blastwaves launched by a strong 
but short explosion (dotted lines). Here the pressure declines 
below P2 after the shock jump, so A//A/ eq rises modestly in- 
ward of w = R s and soon decreases. 

In fact, the right panels of Fig. 1 show the shapes of 
A//A/ eq expected with SZ instruments having two angular 
resolutions; in the top panel we show our results smoothed 
with a Gaussian window of 1" (about 2 kpc at z ~ 0.1), and 
in the bottom panel with 10" (about 20 kpc). The latter reso- 
lution is currently achieved, while the former will be attained 
in the near future, see § 4. 

For frequencies above 220 GHz, similar shapes obtain for 
w < R s , except that the non-thermal contribution is generally 
lower than plotted in Fig. 1. The thermal/non-thermal ratio is 
maximized at v w 370 GHz, to decline only at much larger v. 
At v k> 220 GHz only the non-thermal component survives; 
here the spatial shape of the SZ signal differs significantly on 
two accounts. First, the rise begins at the piston position w = 
R p with a sharp transition of p m \ up to P p ; second, saturation 
occurs at low levels. 

Finally, we discuss why the initial isothermal atmosphere, 
where n oc r~ u applies with w = 2, provides reliable evalua- 
tions of the SZ signals from driven blasts. Flatter runs with 
w < 2 occur in a cluster at r < 100 kpc; here the gravity 
is weaker but the mass m(< r) swept up by the blast in- 
creases faster, and in the shell approximation Eq. (2) yields for 
< P > /P2 a value 6/7 of the isothermal case. Steeper atmo- 
spheres corresponding to 2 < u < 2.5 occur for r > 500 kpc; 



these also yield self-similar blasts propagating at constant M. 
when acted upon by a fading AGN with integrated energy 
AE oc f 2 * 5 " 2 ")/". The outcomes include larger pressure jumps 
but also flatter runs of P(r) at equal values of AE/E (see Lapi 
et al. 2005); overall, somewhat higher levels of P p obtain. 

4. DISCUSSION AND CONCLUSIONS 

SZ signals are widely measured in clusters at levels y w 
10" 4 , see Reese et al. (2002) and Birkinshaw (2004); Lapi, 
Cavaliere & De Zotti (2003) discuss sub-arcmin resolutions 
to detect the integrated effects of AGN outbursts on the ICM 
of clusters and groups. Here we have focused on SZ signals 
resolved at levels of 10" or better, to probe the structure and 
the dynamics of the shocks and blasts so produced. 

The SZ probe is best used in scanning the ICM around den- 
sity jumps selected in X-rays. In fact, the bremsstrahlung sur- 
face brightness proportional to n 2 is well suited for pinpoint- 
ing density jumps and providing positions and Mach num- 
bers of candidate shocks. But measuring T from X-ray spec- 
troscopy must contend with paucity of photons and the narrow 
post-shock range where T exceeds the unperturbed value. 

So in conditions of low surface brightness (outskirts or dis- 
tant structures) the SZ effect will lend a strong hand by un- 
veiling the other key observable, namely, the overpressures 
behind the shocks; this is due to three circumstances. First, 
pressures are sensed directly through the parameter y oc P. 
Second, y is independent of z for sources wider than the in- 
strumental beamwidth. Third, we expect (see Fig. 1) thermal 
pressures to rise throughout a blast continuously driven over a 
crossing time; from shock to piston at radii 40-30% smaller, 
P(r) rises up to values P p w 8 — 9 Pi, considerably larger than 
the shock j ump P 2 /Pi w 2.6-4.8 for M w 1.5-2; this is due 
to the dynamical stress in running blasts. 

This rising behavior of P(r) is just opposite to the run down 
from shock to center expected for free blasts launched by 
short-lived AGN activity. We compare in Fig. 1 the results 
we expect; they constitute an aim particularly interesting for 



4 



A. Cavaliere & A. Lapi 



shocks of intermediate M. pinpointed in close proximity to a 
radiovolume or around a currently shining AGN. Clear study 
cases will be provided by clusters or groups in quiet condi- 
tions, with no sign of outer merger induced dynamics. 

What is needed for SZ probing at z ~ 0.1 is a resolution 
around 10" in the upper /iwave band, already approached with 
the Nobeyama radiotelescope (Kitayama et al. 2004, also 
|http : / /www . nro . nao . ac . jp/ index- e . html/| l. 
Upcoming instruments such as CARJV1A 
Qhttp : / /www . mmarray . orq/| l will do better; at 
resolutions of a few arcsecs, blasts in clusters at z ~ 0.5 
may be probed at levels comparable to the top right panel 
of Fig. 1. ALMA Qhttp : / /www . alma . nrao . edu7) 
with its planned resolutions around 1" and sensitivities down 
to 1 fiK, will do better yet both in the /zwave and in the 
submm band. Following the noise analysis by Pfrommer et 
al. (2005), a 5a detection with ALMA of the thermal SZ 
signal we focus on will require scanning a limited area in the 
vicinity of an X-ray preselected position for a few hours per 
cluster. 

We add that once kpc scales will be resolved, two-fluid 
effects will be interesting as may arise from disequilibrium 
between the electrons and ions (see Ettori & Fabian 1998; 
Fabian et al. 2006). They may cause a gradual rise of the 
electron pressure in front of the shock, to converge behind it 
with the declining ion pressure toward the equilibrium values 
considered above. 

Here we stress three issues. First, how common are shocks 
in clusters? The frequent occurrence of quenched cooling 
flows argues for widespread shocks driven over some Gyrs by 
central AGNs injecting energies AE in excess of the cooling 
losses (see Birzan et al. 2004; Nulsen et al. 2005a, b). 



Second, in what prevailing mode does an AGN release the 
driving energy AE, in mechanical form and radio jets or in 
radiation-driven superwinds? To a first approximation the 
mode little affects the total energy injected, based on the sim- 
ple rule AE ~ afr]c 2 M,At ~ const that we extract from 
Churazov et al. (2005); that is, mechanical energy and jets 
gain on grounds of coupling efficiency (af ?s 1/2 vs. a few 
%) what they lose to radiation on grounds of rjM,. Con- 
versely, the shock energetics alone will not distinguish the 
mode. 

Third, what other probe may help? The specific SZ probing 
we propose can trace the distinguishing features of the injec- 
tion: its timescale (past or ongoing), and its prevailing content 
(relativistic particles or photons). The SZ effect resolved at 
levels of 10" or better will directly detect the mark of a blast 
launched or driven by a powerful AGN, namely, the hydro 
overpressure jumping up at the shock and sustained through- 
out the blast. But when the radiosource drive persists over the 
transit time the pressure actually rises in the blast, and the SZ 
signal is enhanced; it is boosted up by plasma in the cavity if 
the drive is helped by AGN superwinds. 

These outcomes are independent of model details, rather 
they depend on a few overall parameters: the relative injec- 
tion energy AE/E, evaluated from M.; the active time of 
the AGN compared with the blast crossing time; the injec- 
tion mode, whether dominantly in mechanical or radiative en- 
ergy. We conclude that SZ measurements concurring with the 
X-ray imaging can effectively probe the injection mode, the 
dynamics of the blasts, and the history of the driving AGN 
sources. 

We thank our referee for helpful comments. 



REFERENCES 



Birzan, L., Rafferty, D.A., McNamara, B.R., Wise, M.W., & Nulsen, P.E.J., 

2004, ApJ, 607, 800 
Birkinshaw, M., 2004, Carnegie Observatories Astrophysics Series Vol. 3 

(Cambridge: Cambridge Univ. Press) 
Cavaliere, A., Lapi, A., & Menci, N., 2002, ApJ, 581, LI 
Cavaliere, A., & Messina, A. 1976, ApJ, 209, 424 

Churazov, E., Sazonov, S., Sunyaev, R., Forman, W., Jones, C, & Bohringer, 

H., 2005, MNRAS, 363, L91 
Colafrancesco, S., 2005, A&A, 435, L9 
Ettori, S., & Fabian, A.C., 1998, MNRAS, 293, L33 
Fabian A.C., et al., 2006, MNRAS, 366, 417 
Ferrarese, L., 2002, ApJ, 578, 90 
Forman, W., et al. 2005, ApJ, 635, 894 

Heinz, S., Reynolds, C.S., & Begelman, M.C., 1998, ApJ, 501, 126 
Kitayama, T., et al., 2004, PASJ, 56, 17 

Lamers, H.J.G.L.M., & Cassinelli, J.P., 1999, Introduction to Stellar Winds 

(Cambridge: Cambridge University Press), Ch. 12 
Lapi, A., Cavaliere A., & Menci, N., 2005, ApJ, 619, 60 
Lapi, A., Cavaliere, A., & De Zotti, G. 2003, ApJ, 597, L93 



Mazzotta, P., Branetti, G., Giacintucci, S., Venturi, T., & Bardelli, S., 2004, 

Jour, of the Korean Astron. Soc., 37, 381 
McNamara, B.R., Nulsen, P.E.J., Wise, M.W., Rafferty, D. A., Carilli, C, 

Sarazin, C.L., & Blanton, E.L., 2005, Nature, 433, 45 
Nulsen, P.E.J., Hambrick, D.C., McNamara, B.R., Rafferty, D., Birzan, L., 

Wise, M.W., & David, L.P., 2005a, ApJ, 625, L9 
Nulsen, P.E.J., McNamara, B.R., Wise, M.W., & David, L.P., 2005b, ApJ, 

628, 629 

Ostriker, J.P, & McKee, C. 1988, Rev. Mod. Phys., 60, 1 

Parker, E.N. 1963, Interplanetary Dynamical Processes (New York: Wiley) 

Pfrommer, C, EnBlin, T.A., & Sarazin, C.L., 2005, A&A, 430, 799 

Rephaeli, Y, 1995, ARAA, 33, 541 

Reese, E.D., et al., 2002, ApJ, 581, 53 

Scheuer, P.A.G, 1974, MNRAS, 166, 513 

Sedov, L.I., 1959, Similarity and Dimensional Methods in Mechanics (New 

York: Academic Press) 
Stockton, A., Fu, H., Henry, J.P, & Canalizo, G., 2006, ApJ, 638, 635 
Sunyaev, R.A., & Zel'dovich, Ya.B., 1972, Comm. Astroph. Sp. Sc., 4, 173 
Tremaine, S., et al. 2002, ApJ, 574, 740 



